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We investigate some properties of fiat cosmological mod- 
els with a A term that decreases with time as A a a~ m 
(a is the scale factor and m is a parameter < m < 3). 
The models are equivalent to standard cosmology with mat- 
ter and radiation plus an exotic fluid with the equation of 
state p x = (m/3 — l)Px- We study the effect of the decay- 
ing A term on the cosmic microwave background anisotropy 
and by using a seminumeric method we compute the angu- 
lar power spectrum (up to I = 20) for different values of m 
and QmO- We also investigate the constraints imposed on the 
models by the magnitude-redshift test in which high-redshift 
type la supernovae (SNe la) are used as standard candles. 
We obtain the 95.4%, 90%, and 68% confidence levels on the 
parameters m and Q m o and compare them with those arising 
from lensing statistics. Our analysis reveals that the SNe fa 
constraints are stronger for low values of m and fl m o, while 
those from lensing statistics are more important for m ~ f . 
Models with Q m o ~ 0.2 and m ~ f .6 are in good agreement 
with the data. 

PACS number(s): 98.80.Hw 



I. INTRODUCTION 

Flat cosmological models with a cosmological constant 
are currently serious candidates to describe the dynam- 
ics of the universe. There are three main reasons for the 
present interest on these models. First, they can reconcile 
inflation with dynamic estimates for the density parame- 
ter (f2). Observations indicate f2 m o = 0.2 — 0.4 for matter 
that clumps on scales 20 — 30/i _1 Mpc, while inflationary 
models usually predict Qtotal = 1. Second, these models, 
when normalized by data from the Cosmic Background 
Explorer (COBE) predict less power in the perturbation 
spectrum at small scales than standard cold dark matter 
(CDM), in accordance with observations |Q. The third 
motivation for introducing a cosmological constant is the 
"age crisis." In flat models with f2 m o = 1 only if h < 0.57 
(h is the present value of the Hubble parameter in units 
of 100 km/s Mpc -1 ) is it possible to get theoretical ages 
for the universe that are higher than the lowest values 
(~ 12 Gyear) estimated for the galactic globular cluster 
system H]. However, current estimates based on obser- 
vations of Cepheids stars in the Virgo cluster and type 
la supernovae indicate higher values for h || . 

In the past the cosmological constant was introduced 
and, with the improvement of the observational data, 



later discarded. Now, however, the situation may change. 
By taking into account the vacuum contribution to the 
energy-momentum tensor we can define an effective cos- 
mological constant (A e //) that is the sum of two terms, 
the bare cosmological term and 8TvGp vac (Pvac is the vac- 
uum energy density). From quantum field theory we 
should expect p vac ~ ^pi (Mpi is the Planck mass), or 
perhaps another energy density related to some sponta- 
neous symmetry-breaking scale such as Msusy or M wea k 
to the fourth power. The problem lies in that these values 
are enormous when compared with astronomical bounds 
for PA ef f ■ Extreme fine-tuning between &irGp va c and the 
bare A is necessary to make theory compatible with ob- 
servations. 

One possible explanation for a small A term is to as- 
sume that it is dynamically evolving and not constant, 
that is, as the universe evolves from an earlier hotter and 
denser epoch, the effective cosmological term also evolves 
and decreases to its present value [Q. There are also 
strong observational motivations for considering cosmo- 
logical models with a decreasing A term instead of a con- 
stant one. Usually in a dynamical-A cosmological model, 
the distance to an object with redshift z is smaller than 
the distance to the same object in a constant-A model 
with the same value of the density parameter. As we shall 
discuss in Sec. IV, this implies that constraints coming 
from lensing statistics and from high redshift supernovae 
can be considerably weaker in these models jl8| . 

Recently § we suggested a class of models in which 
A decreases as A oc a - '™ [here a is the scale factor of 
the Friedman-Robertson- Walker (FRW) metric and m 
is a constant (0 < m < 3)]. Although Chen and Wu 
H gave some interesting arguments favoring the special 
value to = 2, it is clear that the above functional depen- 
dence with the scale factor is only phenomenological and 
does not come from particle physics first principles. How- 
ever, we believe it deserves further investigation for the 
following reasons. First, these models generalize several 
other models present in the literature. So, by investigat- 
ing their properties we are studying at once the models 
they generalize. Second, since to < 3 the universe age 
in these models is always larger than the age obtained in 
the standard Einstein-de Sitter cosmology and if to < 2 
the age is larger than the one we get in an open model 
with the same f2 m o- This is an important aspect if we are 
interested in solving the "age problem." Further, since 
to < 4 the A term is generally not important during the 
radiation-dominated phase and nucleosynthesis proceed 
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as in the standard model TTj. Finally, the models are 
mathematically simple and in most cases can be treated 
analytically. 

This paper is organized as follows. In Sec. II the ba- 
sic equations of the models are obtained and our main 
assumptions presented. In Sec. Ill we investigate the ef- 
fect of the A term on the cosmic microwave background 
(CMB) anisotropy and compute the angular power spec- 
trum for small values of I for different values of the pa- 
rameter m and Q m o- Constraints on the models from 
high redshift SNe la and from lensing statistics are ob- 
tained in Sec. IV. 



II. THE MODELS AND THE FIELD EQUATIONS 



The above equations are quite general and apply for a 
broad spectrum of models. For instance, if m = 0, the 
usual flat FRW model with a cosmological constant is 
obtained. If we take m = 2, Eqs. (2.3) and (2.4) assume 
the same form of the Einstein equations for open models 
and also appear in some string-dominated cosmologies 
Further, we would obtain the same Einstein equa- 
tions if, instead of a cosmological term, we would have 
considered (beside conserved matter and radiation) an 
exotic x fluid with equation of state, p x = — l) p x 
H . All we discuss here also applies for these cosmologies 
and to emphasize this point in the following we shall use 
the expression "x component" to interchangeably desig- 
nate the A term or the x fluid. 



In this paper we consider spatially flat, homogeneous, 
and isotropic cosmologies with a time-dependent A term: 

A = 8nGp v = 3aa~ m . (2.1) 

Here the parameters a and m are restricted to the range 
0<m<3,a>0 and the factor 3 was only introduced 
for mathematical convenience. We consider the cosmic 
fluid to be a mixture of nonrelativistic matter and radi- 
ation (p r — \p r ) with a perfect fluid energy momentum 
tensor, 

T\ = T/ v + T m % = diag (p, -p, -p, -p), (2.2) 

where p = p r + p m is the total energy density (radiation 
plus nonrelativistic matter) and p — p r is the thermody- 
namic pressure. 

As in Ref. Q we assume that vacuum decays only into 
relativistic particles, such that the nonrelativistic matter 
energy momentum tensor is conserved (p m cx a~ 3 ). The 
radiation energy density has two parts: one conserved, 
^foHq (a /a) 4 , Q r0 = 4.3 x 10~ 5 /i~ 2 , and a second one, 
87rJ(4-m) a ~ m ' wn i cn arises due to the vacuum decay. 
Here oo is the present value of the scale factor and Hq 
is the present value of the Hubble parameter. In the 
following, subscripts will always indicate present values. 

The Einstein equations for the models we are consid- 
ering reduce to two equations: namely, 

GO = °™° ff o(^) +fWf 2 (^) 

+ n xQ H 2 Q ^) m (2.3) 

and 

- = — TT^mO-f^o f ~1 ~ ^toHq ( — ) 
a 2 \ a / Vet/ 

+ ^n x0 HZ^) m , (2.4) 
where VL mo is the matter density parameter and Q x o = 

AaH- 2 g- m 
(4-m) ' 



III. THE ANGULAR POWER SPECTRUM 

The evolution of perturbations in the models we are 
considering may be studied considering two different 
phases ||. For a < clm — ao(j^) rr ™', the energy den- 
sity of the universe is dominated by radiation and/or 
matter, and the contribution from the x component may 
be neglected. During the second phase, characterized 
by a > a,M , the energy density of the universe is domi- 
nated by nonrelativistic matter and/or the x component. 
In this phase, the x component contribution becomes 
more and more important and the deviations from the 
usual matter-dominated cosmological model increase as 
the universe expands. 

To describe the perturbation growth in these two 
phases, we consider the usual linear perturbation the- 
ory. Each Fourier component of the perturbation, 5k (t) , 
grows independently of the other modes and may be re- 
lated to the primordial power spectrum with the help of 
the transfer function: 

S k (t)=T(k)5 k (U), (3.1) 

where 5k(U) is the primordial spectrum, usually taken to 
be 

\5 k {U)\ 2 =Ak n . (3.2) 

In (3.1) T(k) is the transfer function, which incorpo- 
rates all deformations undergone by the fc-mode pertur- 
bation and t may be any time t > U. If the initial pertur- 
bation is adiabatic and n = 1 we have a primordial scale 
invariant adiabatic perturbation that is usually predicted 
by inflationary models. In fact, for flat models without a 
cosmological constant it can be shown that, in this case, 
the quantity k 3 \ 5k | 2 at t — t enter , is independent of k 
for any k. Here, t en t er is the instant when the pertur- 
bation crosses the Hubble radius. However, this is not 
necessarily true if we have a cosmological term. In || we 
adopted a true scale invariant Harrison-Zeldovich spec- 
trum for the models we are analyzing. Here we shall use 
the primordial scale invariant spectrum as defined above. 
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The power spectrum is defined in the standard way by 
P(k) = | 6k(to) 2 where to is the present time. For flat 
models with A = 0, and after radiation-matter equality 
(a > a eq ), all modes, inside and outside the Hubble ra- 
dius, evolve in the very same way. So, 5k(t), in fact, does 
not experience any deformation for t > t eq and P(k) re- 
flects the power spectrum for any t > t eq . Since in our 
model a eq < ajvf, the evolution during the first phase, 
dominated by radiation and matter, may be described by 
simply taking the standard power spectrum at a = clm, 
with the transfer function as computed by Bond and Ef- 
stathiou [pi: 



T(k) = (i + (afc + (6fc) 3 / 2 + {ckfyy 



(3.3) 



where a = QA^lh 2 )' 1 Mpc, b = 3.0(fi/i 2 )- 1 Mpc, c = 
ll^h 2 )- 1 Mpc, and v = 1.13. 

After radiation and x component equality, a > clm, 
the increasingly importance of the x component energy 
density must be taken into account. To describe the per- 
turbations after an, we consider the evolution of matter 
perturbations in the background dominated by matter 
and the x component. We shall use the approximation 
that the x component is smooth. In this case it can be 
shown [[llj that all modes, inside and outside the Hub- 
ble radius, grow at the same rate as in a flat universe 
dominated by nonrelativistic particles. So again no ex- 
tra deformation is introduced in the power spectrum by 
the presence of the x component and, to study the mul- 
tipole expansion of the power spectrum, we are allowed 
to use the transfer function given by Eq. (3.3). 

Since one of our goals is to understand the effect of the 
decaying cosmological term (x component) on the CMB 
anisotropy, we consider the relation between the mass 
density perturbation and the fractional perturbation to 
the CMB temperature %E-. Almost 30 years ago, Sachs 
and Wolfe |Q obtained the expression relating fluctua- 
tions in the gravitational potential on the last scattering 
surface with CMB temperature anisotropies on large an- 
gular scales (9 ^> 1 deg). For a flat universe their formula 
can be written as [1131 
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(3.4) 



where D(t) describes the time dependence of the growing 
mode of 5 m , S m (x,t) = A(x)D(t), j a is the unit vector 
pointing along a null ray from the observer to the source, 
and where k(x) is given by 



S m (x',t) d 3 x' 
D{t) \x-x'\ 



(3.5) 



The first term in the right-hand side (RHS) of (3.4) 
describes the anisotropy caused by the relative motion 



source observer. The effect of the observer's motion may 
be systematically removed from the experimental data, 
and will not be considered any further. The effect of 
the source motion is a Doppler contribution related to 
velocity perturbations on the last scattering surface and 
cannot be eliminated. The second term in Eq. (3.4) in- 
cludes the usual Sachs- Wolfe effect (SW) and a constant 
part evaluated at the present epoch which does not con- 
tribute to the observed anisotropy. Finally, there is the 
last term, usually called the integrated Sachs- Wolfe ef- 
fect (ISW), which represents the redshift or blucshift of 
the photon energy caused by its traveling through regions 
of space with time-varying gravitational potential. Note 
that during the radiation-matter domination phase the 
growing modes behave as §| 



D{a) = D dec (l 



3adec + 
2ap_a 



3a 
2a eQ 



)• 



(3.6) 



where a eq — ao jp 2 - is the scale factor at matter-radiation 
energy density equality. So, since in flat models, during 
the radiation domination, we have a oc t 1 ' 2 , and during 
matter domination the scale factor grows as a oc t 2 ^ 3 , it 
is easy to see that the contribution of the last term in 
Eq. (3.4) vanishes in both cases. This reflects the fact 
that, during these eras, the gravitational potential is time 
independent. However, during the x component domina- 
tion, the gravitational potential will not be constant and 
the ISW term will no longer be zero. In fact its contri- 
bution may be of the same order of the usual SW term 
for the lower modes in the harmonic expansion. In order 
to study the CMB anisotropies in the models described 
above, we keep three contributions in the RHS of (3.4): 
SW, ISW, and Doppler. 

The full temperature correlation function is defined by 

AT AT 

C{a) = (— (") — (™))rVm=cos(a) (3.7) 

where ( . . . ) denotes an average over all positions x and 
all directions h, rh separated by an angle a. 

To compute the average value in Eq. (3.7), we first 
analyze the ISW contribution |l4j . We have 



ao rao 



C lbw {a) = da da'G(a)G(a') 



d 3 k 
(2^)3 



d 3 k' 5k5* k , , j(k'-k).x \ 

(27T) 3 k 2 k' 2 { 

X (e~ *' c -'"'^ 2 ( a ' a o)+^ / -^^2(a / ,ao)^ „ 



n.m— cos 



(a)- 



Here 



h(a, ao) 



dx 



Ho J a y/flmOX + (1 - Q^ojl 



(3.8) 



=, (3.9) 



and 
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G (a) = -(^-iJ'(3+?(^ 
\ a 2 a 



^ ff 2 fi m oa 2 (-^) 3 , 
(3.10) 



where a d = ao(§^) 1/(3 ~ m) and D' = dD/da. Note that 
if m = and flaowe obtain G(a) = 0. Equation (3.8) 
was obtained by making use of the field Eqs. (2.3) and 
(2.4), with Q r o = 0, and the time evolution differential 
equation for D, 



D + 2-D - knGp m D = 0. 
a 



(3.11) 



By using that, 



( e i(fc'-fc).* ^ o = ±_ j d 3 Xoe i& i .,■■„ 



V 



8\k' -k) 



(3.12) 



and 



-ifc.n/2(a,ao)+ife .rali (a ,ao) 



= Jo 



(/ 2 (a, «o) 2 +/2(a',flo) 



)n.m— cos (a) 
2 



1/2 

- 2/ 2 (a, ao)/ 2 (a',ao)cos(a)) 

where jo(x) = sm ^ , we reduce (3.8) to 

dfc P(fc) '" ao z" 10 



(3.13) 



C 



x jo 



r-ri^w f° da r° da ' G ^ G ^ 

JO z7r K Ja M Ja M 

k ^I 2 (a,a ) 2 + I 2 (a',a ) 2 



- 2 7 2 (a, ao)/ 2 (a', a ) cos(a)) 1/2 , 
where P(fc) = AkT(k) 2 . Analogously we obtain 



(3.14) 



3Hq QmoDdec^Q 



,4o e9 (l + fadec/de,),' 

cZfc P(A0 

^2 "p - - 70 ! 2 ^ ( J ( adec ' a o)) sin(a/2)], 



(3.15) 



where I{a dec , a ) = h(a dec , a M ) + h{dM, ao), with 

dx 



1 f a 
h(a dec ,a) = — I 

H o Ja dec 



^rimox + (i - ^ m0 ) 

The Doppler term is given by 

dk P(k) 



(3.16) 



f 

C D °P(a) = a dec 2 Dj ec / 
Jo 



d_d_ 

dx dy 



2tt 2 k 2 
(jo [A: V^ 2 + V 2 - 2x V cos(a)] ) 



x=y=I(a dBC ,a ) 

(3.17) 



The cross terms reduce to 



C ltsw -» w {a) 
'o° 



2>H^Vl. m oD dec a 3 



.3 /-ao fC 

2-) / daG(a) / 

a eg J a M J 



2a e<z (l + |a dec / 



x jo[k ll(a dec ,ao) 2 + h{a, aof 



dk P(k) 
2^2 fc 2 



1 /2 

2/(a a - ec ,ao)l 2 (a,ao)cos(a)) ] (3.18) 
^Dop-siV/^-j _ 1 3HQil m oD dec D dec aQ 



(«) = ( 

d/c P(fc) d 
2^ k 2 dx 



U dec 2 J 



(jo [A: (x 2 + /(a dec ,a ) 2 



1 /2 \ 

- 2x/(a dec ,a )cos(a)) ] (3.19) 

/ x=I(a dcc ,a ) 

C Dop - ISW (a) = -2a dec D dec H daG(a) 

JaM 

dk P(k) d 



„2 , t(„ „ \2 



27T 2 fc 2 dx 



(jo[k (x 2 + I (a, aof 



2xl(a, ao) cos(a)) 1 ^ 2 ] ) 

/ X = I( 



(3.20) 



x=I(a,dec,ao) 

The total contribution to (3.7) may be written as 
C(a) = C ISW + C sw + C Dop + C ISW - SW 

_|_ QDop—SW _|_ QDop—ISW ^ 22) 

It is convenient to expand the temperature correlation 
function in angular multipoles, using Lcgcndrc Polyno- 
mials p(cos(a)): 



00 9/ 4- 1 



1=0 



(3.22) 



and Ci = 2tt j\dcos(6)Pi(cos(6))C(6). With the help 

of the relation jo[ky/r 2 + q 2 — 2rocos(0)] = X)^o(^ n + 
l)P n (cos(9))j n [kr]j n [kq], where ji[x] is a spherical Bessel 
function, the integrals in the contributions for C; decou- 
ple and the angular power spectrum becomes 



ISW 



with 



where 



dfc 



T ,/ ^ 



P(fc)(A^(fc)) 2 



(3.23) 
(3.24) 



/■ao 

X;(fc)= / doG(a)ji[fc/(a,ao)], (3.25) 

J «M 



SW _ ^ ( 3Ho^rnoD d ec(lo 



tt \4a e9 (l + ^a dec /a eq ) 



f°° dk 

Jo ^ P ik)3i[kl{a d ec,a )l (3.26) 
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x=I(a dec ,a a ) 



(3.27) 



c, 



ISW-SW 



ira eq (l + ^adec/a-eq) 



x J ^P{k)ji[kI{a dec ,a Q )]Xi(k), (3.28) 



x ii[fc/(a dec ,a )]^: (7l[*a;])x=J(o„. e , 



and 



AddecDa 



,ao) 



(3.29) 



-p ( fc)x ; (fc)^^[M)^ (adcc , ao) . 



(3.30) 



Using the expressions above, we compute the angular 
power spectrum for different values of m and il m o- Cal- 
culations are done in a semianalytical way through the 
following steps. To make the integration in Xi(k) [Eq. 
(3.25)], we change variables. By defining w = I (a, do), 
and by using an interpolated polynomial, we rewrite the 
integral in the new w variable as 



X, 



f' 

Jo 



dw CI{w) ji[kw]. 



(3.31) 



Here CI(w) is an interpolated polynomial, which re- 
places the function C(w) = (w))G(I^ (to)), 
where I^~ 1 \w) stands for the inverse of w = I(a,ao). 
The precision of this interpolation may be estimated by 
comparing CI(w) and the original function C{w). In 
the models we consider, this comparison shows that us- 
ing the seventh-order interpolated polynomial, CI(w) — 
12n = o a nW n 7 introduces an error always smaller than 
0.07%. With the help of this approximation, the result- 
ing seven integrals in Xi can be analytically computed. 
The form of these expressions remains the same for all 
models that we want to consider. Changing among mod- 
els only changes the numerical values of a n together with 
the upper limit in the integral, %. Once we have cal- 
culated these auxiliary functions Xi(k), we may compute 
the coefficients in Eqs. (3.24), (3.28), and (3.30) by di- 
rect integration in k. In the limit k — > 0, the integrand 
is well behaved. For large values of k, the integrand de- 
creases very fast. However, while decreasing, these in- 
tegrands oscillate wildly, making the integral procedure 
slower and slower. To avoid this technical difficulty, we 
truncate the integral at some value of k, and we certify 



ourselves that the truncation is affecting our results in a 
controlled way. For values of k larger than the truncation 
value, we compute the integrals replacing the oscillating 
ji(kl) by an envelope function, and the result of these 
integrals is taken as an upper bound of the error intro- 
duced by the truncation in k. This error grows steadily 
with I, starting with 0.03% for I = 2 and growing up to 
the range 2 — 3% for / = 20 and different values of m. 

The numerical results obtained with the method de- 
scribed above are plotted in Figs. 1 and 2. 




0,n m0 = 0.2 
m=1,n m0 = 0.2 
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FIG. 1. The angular power spectrum 1(1 + 1)C\/6C2 x I is 
shown for different values of m and fl m o- 





log 1 

FIG. 2. The Doppler, the Sachs-Wolfe, and the integrated 
Sachs- Wolfe contributions to the angular power spectrum are 
shown. For the figure we considered the special model for 
which m — 2 and f2 m0 = 0.2. 

In Fig. 2 we plot for m — 2 and f2 m o = 0.2 the separate 
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contributions from the Sachs- Wolfe, Doppler and inte- 
grated Sachs- Wolfe effects. Crossed terms are not shown. 
Basically, in the range 2 < I < 20, the Sachs- Wolfe contri- 
bution remains always important, with a slight decrease. 
The integrated Sachs- Wolfe is large for small values of I, 
but soon it becomes unimportant, reducing to just 2% 
of Ci at I = 20. On the other hand, the Doppler con- 
tribution starts very small and grows fast to become the 
dominant contribution at I = 20. 

For the analyzed models, the curves 1(1 + 1)C;/6C2 x I 
in Fig. 1, show a common feature, a minimum for small 
values of I. We see that, for fixed Sl m o and h, the posi- 
tion of this minimum changes with to, becoming deeper 
as to grows larger. This feature reflects the fact that the 
contribution from the ISW effect is larger for larger val- 
ues of the parameter m. We also observe that for fixed 
to and h and different values of the matter content, the 
minimum grows deeper for smaller values of Q m Q. Again 
this feature agrees with the fact that it is the x com- 
ponent term which brings up the nonzero ISW effect, 
which, by its turn, makes the minimum deeper. So, since 
Qtotai = 1, less nonrelativistic matter means more x com- 
ponent contribution, more ISW effect, and consequently 
deeper minimum. 

IV. CONSTRAINTS FROM HIGH-REDSHIFT 
TYPE IA SUPERNOVAE AND LENSING 
STATISTICS 

The Supernovae Cosmology Project is an ongoing pro- 
gram to systematically search and study high-z super- 
novae. In a recent report |lq ] Perlmutter et al. analyzed 
seven SNe (with redshift z — 0.35 — 0.46), of more than 
28 supernovae discovered, and obtained constraints on 
cosmological parameters, specially on the cosmological 
constant. Their preliminary result, Q\ < 0.51 at the 
95% confidence level, strongly constraints models for the 
universe with a cosmological constant. In this section we 
use their observational results and adapt their procedure 
to constraint the class of models described in Sec. II. 

The basic idea is to use type la supernovae as standard 
candles for the classic magnitude-redshift test. As in 
we express the apparent bolometric magnitude m(z) as 

m(z) = M + 5 log di(z,O, m0 ,m) (4.1) 

where in our case the luminosity distance (in units of 
Hq 1 ) is given by, 

di(z, Q m0 ,m) = c(l + z) 

r dz , A . 

Jo y/(i + z)m m0 + (i + z ) m (i - n m0 ) 

In Eq. (4.1) 

M = M-5 logffo + 25 (4.3) 



is the "zero point" magnitude (or Hubble intercept mag- 
nitude) , which is estimated from the apparent magnitude 
and redshift of low-redshift (z < 0.1) SNe la. The nearby 
supernovae data set used by Perlmutter et al. in the de- 
termination of M. were those 18 SNe la, discovered by 
the Calan/Tololo Supernovae Search jl6| for which the 
first observations were made no later than 5 days after 
maximum. 

Although SNe la are very similar explosion events, it 
is now known that they do not constitute a completely 
homogeneous class. Recently progress was made in the 
study of their inhomogeneities. Phillips |17| showed that 
there is a correlation between absolute magnitudes (M) at 
maximum light and the initial decline parameter Atois, 
the B-magnitude decline in the first 15 days after maxi- 
mum. By studying the Calan/Tololo type la supernovae, 
Hamuy et al. |16| confirmed the existence of the correla- 
tion suggested by Phillips and obtained a prescription for 
correcting observed B-magnitudes to make them compa- 
rable to an arbitrary "standard" SNe la light curve of 
width Ami5 = 1.1. After adding the correction to the 
SNe la low-z set they reduced the magnitude dispersion 
from 0.26 to a^ muy = 0.17. Perlmutter et al. did the 

1V1 b ,corr 

same for the seven high-z supernovae and reduced the 
dispersion from 0.27 to 0.19 mag. 




0.25 0.5 0.75 1 1.25 1.5 

m 



FIG. 3. Constraints imposed on the models by the magni- 
tude-redshift test in which high-z SNe la are used as standard 
candles. The 95.4%, 90%, and 68% confidence levels on the 
parameters m and Q m o are shown in the figure. 

In our computations we follow jl5[ ] and use the cor- 
rected B-magnitude intercept at ATO15 = 1.1 mag, 
■M^j 'corr = ™3.32 ± 0.05, and consider only those 5 su- 
pernovae that have ATO45 values in the range 0.8 — 1.5 
mag, which is the range of values investigated in the 18 
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low-z supernovae dataset. To construct the % 2 we used 
the data points outer error bars of Perlmutter et at, 
which are obtained by adding in quadrature the inner 
error bars of niB,corr (the apparent B-magnitude after 
width-luminosity correction) to a^ muy . 

In Fig. 3 we plot the 95.4%, 90%, and 68% confidence 
levels for the parameters to and fl m o- We see that for the 
interesting range fi TO o ~ 0.2, models with to ~ 1.3 are in 
good agreement with the data. The goodness-of-fit (as 
defined in [jDsf) for the considered case is 0.59. 

If we fix m = we recover the result of Perlmutter et 
al. for constant A, fi A = 0.061^34, with ft A < 0.51 at 
the 95% confidence level (one tail). Another interesting 
case is m = 2. In this case we again have to consider 
one degree of freedom in A% 2 = x 2 — Xmin > anc ^ we ^ nc ^ 
!l m o = 0.83+0^9 (l^) w hh goodness-of-fit equal to 0.75. 
So, models with m = 2 and f2 m o > 0.14 are also in good 
agreement with the data. 

It is interesting to compare the above constraints with 
those arising from lensing statistics. In Fig. 4 we plot 
contours of constant likelihood (95.4%, 90%, and 68%) 
for the models, arising from lensing statistics. We used 
data from the HST Snapshot survey, the Crampton sur- 
vey, the Yee survey, the ESO/Liege survey, The HST GO 
observations, the CFA survey, and the NOT survey |jl9| . 
We considered 859 (z > 1) quasars plus 5 lenses and to 
obtain the contours we adapted from jfj| the procedure 
denoted there as "A2" . 




FIG. 4. Contours of constant likelihood (95.4%, 90% and 
68%) arising from lensing statistics are shown. 

By comparing Fig. 3 with Fig. 4 we see that the SNe 
la constraints are stronger for low values of to and f2 m o 
(to < 1 and O m o < 0.4) while the lensing ones are more 
important for larger values of the parameter to (to > 1). 



From Fig. 4 we observe that models with O m o ~ 0.2 and 
to ~ 1.6 are in good agreement with the data. 
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